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THE STRUCTURE OF HALOS IN SELF-INTERACTING COLD DARK 

MATTER MODELS 

P. Colin, 1 V. Avila-Reese, 2 O. Valenzuela, 3 C. Firmani 4 
RESUMEN 

Usando simulaciones de alta resolution, hemos estudiado la estructura y subestructura de halos corrcspondientes 
a galaxias tipo Via Lactca y a cumulos de galaxias para una cosmologia ACDM con parti'culas oscuras auto- 
interactuantes. La section eficaz, ctdm-i en nucstros modelos puede ser constante o proporcional al inverso de la 
velocidad relativa, v 12 . Esta ultima dependencia de o^m con v 12 resuelve el problema del "nucleo empinado" 
de la cosmogonia ACDM y predice una poblacion de subhalos similar a la que se observa en los halos CDM. 

ABSTRACT 

High-resolution numerical simulations were performed to study the structure and substructure of Milky Way- 
and cluster-sized halos in a ACDM cosmology with self-interacting (SI) dark particles, where the particle cross 
section, <JdMi is assumed constant or inversely proportional to the relative velocity. We conclude that the 
cuspy halo problem at galaxy scales of the ACDM cosmogony can be solved in the latter case. At the same 
time, the subhalo population remains roughly similar to that seen on CDM halos. 
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1. INTRODUCTION 

The potential problems of the popular ACDM 
model at small scales have motivated the search for 
modifications that keep the succesful predictions of 
the model at large scales unaltered. Two of these 
problems are that (a) ACDM halos seem to be more 
concentrated and cuspy than what is suggested by 
the rotation curves of dwarf and LSB galaxies (see 
references below) , and that (b) the number of subha- 
los predicted by the ACDM model on galactic scales 
overwhelms the observed one (Klypin et al. 1999). 
While the latter problem has a well-motivated as- 
trophysical solution (see e.g., Benson et al. in this 
volume), the former one seems to require a solution 
which goes beyond the standard model of structure 
formation. A posible solution for the "cuspy" prob- 
lem was proposed by Spergel & Steinhard (2000) (see 
also Firmani et al. 2000) who introduced the con- 
cept of self-interacting (SI) DM particles within this 
context. The SIDM cross section per unit of par- 
ticle mass, ctdm-i is a free parameter which may be 
constant or dependent on the relative velocity of the 
colliding particles, v 12 : <Jdm = a o^ioo~ a ^ where v 10 o 
is V12 in units of 100 kms . 

A number of authors have constrained the range 
of the pair (er , a) of values and concluded that the 
relevant regime for structure formation is the opti- 
cally thin one (small cross sections; see e.g., Hennawi 
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& Ostrikcr 2002 and referenes therein). By means 
of cosmological simulations, we explore in more de- 
tail the effects of varying (<7Q,a) on halo structure 
and subhalo population (Colin et al. 2002, hereafter 
C2002 ). An outline of the results is presented in 
§3. In order to explore the viability of the SIDM 
model, it is important to establish the observational 
constraints. This is done in the next section. 

2. OBSERVATIONAL CONSTRAINTS 

Several recent observational studies, using the 
highest sensitivity and spatial resolution in HI, Ha 
and CO lines, show that halos of dwarf and LSB 
galaxies seem to be less concentrated in the center 
than what is predicted by the ACDM model (e.g., 
Cote, Carignan, & Freeman 2000; Blais-Ouellette, 
Amram, & Carignan 2001; Bolatto et al. 2002; Am- 
ram & Garrido 2002; de Blok et al. 2001a, 2001b; 
Marchesini et al. 2002; see also the contributions by 
Bosma and by de Blok in this volume). The inner 
slopes inferred for the halo density distribution are 
typically around —0.5, shallower than the slopes of 
ACDM halos. There are also pieces of evidence of 
soft cores in normal disk galaxies (e.g., Corsini et 
al. 1999; Firmani & Avila-Reese 2000; Borrielo & 
Salucci 2001; Salucci 2001) and in elliptical galaxies 
(Keeton 2001), although the evidence is less direct 
than in the case of dwarf and LSB galaxies. 

An important question is that of the scaling laws 
of halo cores. Assuming the non-singular isothermal 
or pseudo- isothermal halo models, several authors 
have tried to find this law. For a sample of dwarf and 
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LSB galaxies with high-quality rotation curves and 
two clusters of galaxies studied with gravitational 
lensing, Firmani ct al. (2000, 2001) have found that 
the halo central density, when plotted against the 
maximum circular velocity, v max , exhibits a large 
scatter with no correlation with v max ; the core ra- 
dius, on the other hand, tend to correlate with v max 
with a slope smaller than 1. Other workers have ar- 
rived to similar conclusions (see references in C2002). 
The data on clusters have increased in the last two 
years, mainly from high-resolution Chandra X-ray 
studies. The mass distributions for most of the X- 
ray clusters are well fitted by both the NFW and the 
pseudo-isothermal profiles. The studies with strong 
gravitational lensing favour shallow halos. 

For a sample of dwarf and LSB galaxies with 
high-resolution rotation curves and for X-ray and 
lensing clusters of galaxies, we infer that the cen- 
tral density p c ,-i, defined as the density where 
the logarithmic slope of the density profile be- 
comes lower than —1, has values between ~ 10 16 — 
10 17 h 2 M Q Mpc~ 3 from dwarf- to cluster-sized ha- 
los, with no evidence of a dependence on v max . The 
predictions of the ACDM model for p c -i for cluster- 
sized halos are within this range, but for galaxy-sized 
halos, Pc,-i is larger than 10 17 h 2 M Q Mpc -3 . Can 
the SIDM cosmology solve the discrepancy? 

3. MODELS AND RESULTS 
The Adaptive Refinement Tree (ART) N-body 
code (Kravtsov, Klypin, & Khokhlov 1997) has been 
used to run the N-body simulations. A pair of 
particles collide if the distance d = — Aln(l — P) 
(A = l/ptiDM is the mean free path and P is a ran- 
dom number distributed uniformly between and 1) 
becomes lower than the distance v^Af, where V12 is 
the relative velocity between the particle and one of 
its nearest neighbors, and At is the time step. The 
ACDM model was used throghout. 

The core evolutionary behavior of NFW halos 
with the self-interaction on is shown in our exper- 
iments with monolithic halos. Firstly, the core ex- 
pands due to the heat inflow from the hotter sur- 
roundings (for CDM halos, the 3D velocity disper- 
sion decreases towards the center). There is a radius 
below which quantities like the total energy or the 
heat capacity, C, become positive. Therefore, the 
heat inflow leads to the isothermalization of the core. 
Secondly, after the maximum expansion of the core, 
the region where C > decreases, heat starts flow- 
ing from a zone where C is actually now negative. 
Thirdly, the gravothermal instability is triggered and 
the core collapses. The system moves from a mini- 
mum to a maximum entropy state. The time scales 



of these processes for monolithic halos depend on 
&dm ■ For cosmological halos these time scales can 
be modified by extra dynamical heating (e.g., merg- 
ers): the halo periphery can be heated and the tem- 
perature there can become higher than in the core, 
the heat flows to the core and it expands until the 
overall halo isothermalizes. This process can delay 
the core collapse, but will not reverse the gravother- 
mal instability to a runaway core expansion. 

We have further investigated Milky Way (MW)- 
and cluster-sized halos from cosmological simula- 
tions, using several values for the pair [ao,a). The 
parameter that defines the evolution under SI is the 
number of collisions per particle after a given time, 
N co i oc p <j dm v rms t. From all of our runs, we see 
that if N co i at the present epoch is < 2 — 5, then the 
halo is either in the core-expansion phase or has just 
had the gravothermal catastrophe triggered. The 
halo density profiles are affected (flattened) only in 
the innermost regions, r < 0.05R,ji r . When a dm is 
constant [a. = 0), the cores of low velocity halos are 
on average less influenced by SI than the high veloc- 
ity ones; therefore, p c ,—i will depend on scale (but 
see Dave et al. 2001). When odm <x I/V12 ( a = l)i 
we expect halos of different sizes to have roughly sim- 
ilar core densities (Firmani et al. 2001; Yoshida et 
al. 2001) Nevertheless, these simple reasonings apply 
strictly only for monolithic halos. The cosmological 
mass assembly history may dramatically influence 
the z = structure of halos with a significant cross 
section; for example, for values of oo = 3 cm 2 gr _1 
and a = 0, we find that MW- or cluster-sized halos 
are or are not well inside the core collapse phase by 
2 = 0; this depends essentially on whether halos have 
suffered recent major mergers. We see clearly that 
the smooth mass accretion is not important; only the 
violent mergers can delay the core collapse phase. 

We find that only SIDM models with <tdm oc 
I/V12 are able to predict halo central densities that 
do not depend on scale, as observations suggest. For 
o dm = 0.5 — 1.0(l/vioo)cm 2 gr _1 , dwarf to cluster- 
sized halos have p c ,-\ values within the range in- 
ferred from observations. In Fig. 1, we present 
a comparison between a MW-sized SIDM (a dm = 
0.5(l/vioo)cm 2 gr _1 ) and a CDM halo (the so-called, 
MW-sized, fiducial halo in C2002). Slight differences 
can be noticed from this figure: (a) the core of the 
SIDM halo is less concentrated, the black spot at 
the center of the halo is less dark and appears uni- 
formly spread over a larger area, (b) small subhalos 
are slightly more numerous in the SIDM halo and 
have a different radial distribution, and (c) the SIDM 
halo seems to be rounder (a quantitative and more 



SELF-INTERACTING DARK MATTER 



3 




Fig. 1. Distribution of dark matter particles inside a box of 0.5 h _1 Mpc on a side centered on a MW-sized halo with 
(tdm = (left, CDM) and with ctdm = 0.5(l/vioo) (right, SIDM). Particles are color-coded on a gray scale according 
to their local density. The local density at the particle positions was computed using SMOOTH, of Astronomy of the 
University of Washington. 



detailed comparison can be found in C2002). 

In general, we find that the number of subhalos 
within the SIDM halos simulated here is largely sup- 
prcscd only for (ao,a) = (3.0,0.0). Subhalos may 
survive a longer time in SIDM halos than in CDM 
ones because the inner tidal force in the former is 
weaker than in the latter. The structure and popu- 
lation of subhalos are determined by the interaction 
between the hot host halo particles and the cooler 
subhalo particles, rather than by internal processes 
in the subhalos. Overall, SIDM subhalos become 
puffier than their CDM counterparts; they are com- 
pletely evaporated only in the limit of high a dm- 

Our simulations also show that SI start flattening 
the inner density profiles of growing CDM halos since 
early epochs in such a way that the growth of super 
massive black holes by accretion of SIDM (Hcnnawi 
& Ostriker 2001) does not seem to be an efficient 
process, and therefore can not be used as a criterion 
to constrain SIDM models. 

This work was supported by CONACyT grants 
33776-E to V.A. and 36584-E to P.C. 
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